theta = 0:0.1:4*pi;
y  = sin(theta);
maxttt = max_cruvature(theta,y)

